Classical confinement of test particles in higher-dimensional models: stability criteria 

and a new energy condition 
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We review the circumstances under which test particles can be localized around a spacetime sec- 
tion Eo smoothly contained within a codimension-1 embedding space M. If such a confinement is 
possible, Eo is said to be totally geodesic. Using three different methods, we derive a stability con- 
dition for trapped test particles in terms of intrinsic geometrical quantities on Eo and M; namely, 
confined paths are stable against perturbations if the gravitational stress-energy density on M is 
larger than that on Eo, as measured by an observed travelling along the unperturbed trajectory. We 
confirm our general result explicitly in two different cases: the warped-product metric ansatz for 
(n + l)-dimensional Einstein spaces, and a known solution of the 5-dimensional vacuum field equa- 
tion embedding certain 4-dimensional cosmologies. We conclude by defining a confinement energy 
condition that can be used to classify geometries incorporating totally geodesic submanifolds, such 
as those found in thick braneworld and other 5-dimensional scenarios. 

PACS numbers: 04.20.Jb, ll.10.kk, 98.80.Dr 



I. INTRODUCTION 

The past half-decade has seen a notable upswing in 
interest in non-compact higher-dimensional theories of 
physics. Most of this attention can be attributed to re- 
cent advances in string theory, which have postulated 
that we are living on a (3 + l)-dimensional hypersur- 
face embedded within some higher-dimensional manifold. 
Such "braneworld" scenarios have been extensively ana- 
lyzed in the literature, and have been used to address 
issues such as the hierarchy problem of particle physics 
ELQUiQ) as well as the idea that the post- inflationary 
epoch of our universe was preceded by the collision of 
D3-branes pi El- In all fairness, it should be mentioned 
that the current flurry of interest in braneworld scenarios 
has been preceded by numerous other models making use 
of large or infinite extra dimensions 0, H, El El • 

In some braneworld scenarios, the idea of non-compact 
extra dimensions is made more palatable by postulating 
that the particles and fields of the standard model are 
confined to the brane universe. If we adopt the most 
conservative point of view, the notion of confinement is a 
prerequisite for any reasonable theory with non-compact 
extra dimensions; without such an assumption, the fact 
that we so not commonly see objects flying off in un- 
seen directions becomes a thorny issue. In the context 
of a particular strin g th eory-inspired model put forth 
by Horava & Witten [la . Il4j , lower-dimensional confine- 
ment is a natural consequence of the idea that standard 
model degrees of freedom are associated with open strings 
that have endpoints residing on a Dp-brane. Conversely, 
since gravitational degrees of freedom are associated with 
closed strings, the graviton in such models is assumed 
to propagate both in the bulk and on the brane. Phe- 



nomenological 5-dimensional realizations based this idea 
model the brane as a 4-dimensional domain wall or de- 
fect y, 3- The discontinuity in the 5-geometry about 
the defect forces the graviton ground state to be sharply 
localized on the brane, which allows for the recovery of 
standard Newtonian gravitation in the low-energy limit. 
This kind of localization extends to other types of fields, 
thus representing a sort of concretization of the confine- 
ment mechanism envisioned in the original string model. 
In addition, if the matter localized on the brane satisfies 
the appropriate energy conditions and the Z2 symmetry 
is obeyed, one can show that test particles can be grav- 
itationally confined to a small region about the defect 
[l5|. This acts as a classical confinement mechanism. 

A natural generalization of models involving thin ge- 
ometric defects is scenarios involving thick, smooth do- 
main walls [lii. fl7j |. There are a couple of reasons to 
consider such models: First, there is a natural minimum 
length in string theory given by the string scale, so the 
idea of an infinitely thin geometric feature is somewhat 
suspect even in a phcnomenological model. Second, one 
would like to see these braneworld scenarios resulting 
from some genuine solutions of supergravity, which are a 
priori smooth and differentiable manifolds. The question 
is: what becomes of the confinement paradigm in bulk 
manifolds devoid of thin domain walls? For test particles 
in scenarios with one extra dimension, the answer is well 
known: if the brane has vanishing extrinsic curvature, 
geodesies may be naturally hypersurface-confined with- 
out the invocation of external non-gravitational forces. A 
surface with zero extrinsic curvature is sometimes called 
totally geodesic. 1 But what of the stability of the trajec- 
tories confined on these surfaces? That is, if one perturbs 
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1 An alternative name for a totally geodesic submanifold 
"geodesically complete." 
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a confined trajectory slightly off of a totally geodesic sub- 
manifold, will it naturally return or not? In other words, 
under which conditions is a totally geodesic hypersur- 
face gravitationally attractive? For obvious reasons, such 
questions are of direct relevance to any serious attempt to 
classically describe our universe as a smoothly-embedded 
hypersurface on which we are gravitationally trapped. It 
is possible that this classical stability issue is irrelevant 
at the quantum level — perhaps because stable particle 
confinement can be guaranteed by other means — but 
for the purposes of this study we will assume that the 
classical formalism is applicable. 

In this paper, we propose to address these issues n- 
dimensional totally geodesic submanifolds smoothly em- 
bedded in a space of one higher dimension, with cither 
timelike or spacelike signature. We will utilize quite gen- 
eral methods that will ensure our results apply to any 
geometry and choice of coordinates in the bulk or on 
the submanifold. In Section [Hj we describe our geomet- 
ric construction. In Section IlIIl we review the covariant 
splitting of test particle equations of motion developed in 
(15L llSl Il9| and use it to derive the zero-extrinsic curva- 
ture condition for totally geodesic submanifolds. Then, 
we find the stability condition for the confined trajecto- 
ries, which is that the double contraction of the particle's 
velocity with the Ricci tensor of the bulk is greater than 
the double contraction with the Ricci tensor of the sub- 
manifold. In more physical terms, the stability of trapped 
particles demands that the locally measured gravitational 
density of the bulk is bigger than the density of the effec- 
tive lower-dimensional matter living on the brane. 2 We 
briefly discuss the nature of the latter, emphasizing that 
the stress-energy content of the submanifold — as per- 
ceived by an observer ignorant of an extra dimension — 
is made up from contributions from the "real" higher- 
dimensional matter as well as the bulk Weyl tensor. For 
good measure, we derive the stability condition using two 
additional methods: the geodesic deviation equation in 
Section llVAl and the Raychaudhuri equation in Section 
IIVBI We confirm the correctness of our general result for 
the special case of the warped-product metric ansatz in 
Section IV Al and consider a simple 5-dimensional model 
of the solar neighborhood. In Section IV Bl we show that 
our stability condition is also correct in the Liu-Mashoon- 
Wesson solution |20l l2ll| of the 5-dimensional vacuum 
field equations. Section IVII summarizes our work and 
presents the confinement energy condition, which ensures 
that all timelike trajectories on a totally geodesic sub- 
manifold in a given bulk geometry will be stable. This 



energy condition can be used to classify solutions of the 
thick braneworld on other 5-dimcnsional scenarios. 

Conventions. Uppercase Latin indices run from 
to n, while lowercase Greek indices run from to 
n — 1. Higher-dimensional curvature tensors are dis- 
tinguished from their lower-dimensional counterparts by 
hats. Higher and lower dimensional covariant differentia- 
tion operators are denoted by or V a , respectively. A 
center dot indicates the scalar product between higher- 
dimensional vector fields; i.e., u ■ v 



u a va- 



II. GEOMETRIC CONSTRUCTION 

We will be concerned with an (n+l)-dimensional man- 
ifold (.M, <?ab) on which we place a coordinate system 
x — {x A }. Sometimes, we will refer to M as the "bulk 
manifold." In our working, we will allow for two possibil- 
ities: either there is one timelike and n spacelike direc- 
tions tangent to M, or there are two timelike and (n — 1) 
spacelike directions tangent to M. Hence, the signature 
of gAB is 

sig(<?A B ) = (- + ••• + e), (1) 
where e = ±1. We introduce a scalar function 



t(x), 



(2) 



which defines our foliation of the higher-dimensional 
manifold with the hypersurfaces given by I — constant, 
denoted by S^. If there is only one timelike direction 
tangent to M, we assume that the vector field n A nor- 
mal to is spacelike. If there are two timelike directions, 
we take the unit normal to be timelike. In either case, 
the submanifold tangent to a given hypersurface con- 
tains one timelike and (n — 1) spacelike directions; that 
is, each hypersurface corresponds to an n-dimensional 
Lorentzian spacetime. The normal vector to the T,i slic- 
ing is given by 



tla — 8a£, n ■ n 



(3) 



The scalar $ which normalizes n A is known as the lapse 
function. We define the projection tensor as 



flAB = gAB - EUATlB ■ 



(4) 



This tensor is symmetric {Hab = hsA) and orthogonal 
to ha- We place an n-dimensional coordinate system on 
each of the hypersurfaces y = {y a }- The n holonomic 
basis vectors 



2 Roughly speaking, the gravitational density of a given matter- 
energy distribution differs from the ordinary density by terms 
involving the pressure. For example, according to an observer 
comoving with a (n + l)-dimensional perfect fluid, the gravita- 
tional density — as we define it below — is [(n — 2)p+np]/(n — 1). 
It is important because it appears naturally in the Raychaudhuri 
equation, as we will see in Section llV Bl 



dx A 
dy a 



= 



(5) 



are by definition tangent to the £^ hypersurfaces and or- 
thogonal to n A . It is easy to see that e A behaves as a 
vector under coordinate transformations on M [<fi : x — > 
x(x)\ and a one-form under coordinate transformations 
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on E^ [if) : y — > We can use these basis vectors 

to project higher-dimensional objects onto E^ hypersur- 
faces. For example, for an arbitrary one-form on M we 
have 



(6) 



Here T a is said to be the projection of Ta onto E^. 
Clearly T Q behaves as a scalar under </> and a one-form 
under tp. The induced metric on the E^ hypersurfaces is 
given by 

h a/ 3 = e A e B g A B = e A e B h A B- (7) 

Just like <?ab, the induced metric has an inverse: 

h^hyp = 6 a p. (8) 

The induced metric and its inverse can be used to raise 
and lower the indices of tensors tangent to E^, and change 
the position of the spacetime index of the e A basis vec- 
tors. This implies 



5 a 



0- 



(9) 



Also note that since Hab is entirely orthogonal to n A , we 
can express it as 

hAB = h a(i e a A e P B . (10) 

At this juncture, it is convenient to introduce our defini- 
tion of the extrinsic curvature K a p of the E^ hypersur- 
faces: 

K a p = e A e B V A n B = \e. A e B £ n h A B- (H) 

Note that the extrinsic curvature is symmetric (K a p — 
Kp a ). It may be thought of as the derivative of the in- 
duced metric in the normal direction. This n-tensor will 
appear often in what follows. 

We will also require an expression that relates the 
higher-dimensional covariant derivative of (n+ l)-tensors 
to the lower-dimensional covariant derivative of the cor- 
responding n-tensors. We have that the n-dimensional 
Christoffel symbols are given by 



1 /9 7 



e B e a A V B ef 



(12) 



s 7 e^vse^. 

This allows us to deduce that for one-forms, the following 
relation holds: 



VpT a = e B e A V B (h A u T c ) 



Cr 



(13) 



where \7 B is the covariant derivative on M defined with 
respect to §ab and is the covariant derivative on E^ 
defined with respect to h a p- The generalization to ten- 
sors of higher rank is obvious. It is not difficult to confirm 
that this definition of V Q satisfies all the usual require- 
ments imposed on the covariant derivative operator. 

Finally, we note that {y, £} defines an alternative co- 
ordinate system to x on M. The appropriate diffcomor- 
phism is 



dx 1 



e A dy a + £ A d£, 



(14) 



where 



£ A = 



( dx J 



V d£ 



(15) 



y a — const. 



is the vector tangent to lines of constant y a . We can 
always decompose higher dimensional vectors into the 
sum of a part tangent to E^ and a part normal to S^. 
For £ A we write 



N a e A + $n J 



(16) 



This is consistent with £ a 8a£ = 1, which is required by 
the definition of £ A , and the definition of n A . The n- 
vector N a is the shift vector, which describes how the y a 
coordinate system changes as one moves from a given E^ 
hypersurface to another. Using our formulae for dx A and 
we can write the higher dimensional line element as 



dsf M) = g A B dx A dx B 

= K p {dy a + N a d£)(dy 13 + N p d£) 
+e<f> 2 d£ 2 . 



(17) 



This reduces to ds 2 Si) — h a pdy a dyP if d£ — 0. It is also 
possible to express the extrinsic curvature in terms of $ 
and N a : 



K, 



a/3 



2$ 



(18) 



where £n is the Lie derivative in the direction of the 
shift vector. 

In this paper, we will be primarily concerned with the 
Gaussian-normal coordinate gauge that has been termed 
canonical by some authors . This is defined by the 
following choices of foliation parameters: 



$ = 1, N a = 0. 



(19) 



Obviously, this choice will result in significant simplifica- 
tion of many of the preceding and following formulae. 



III. CONFINEMENT OF TEST PARTICLES 

The equations of motion for a test particle travelling 
through M are taken to be 



u A V A u B = T l 



u • u = ft, 



dx A 



(20) 



where k = —1, 0, +1 to allow for massive, null and tachy- 
onic particles respectively, A is an affine parameter, and 
T is some non-gravitational force per unit mass. One can 
decompose these equations into relations involving the 
particle's velocity tangent to the E^ foliation u a = e a ■ u 
and parallel to the normal direction u n — n ■ u. This 
was first done in |l8j for a 5-dimensional model with 
a spacelike extra dimension and pure geodesic motion, 
then generalized to accelerated trajectories and an ex- 
tra dimension of arbitrary signature in |l5j| , and further 
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adapted to arbitrary dimension and refined notation in 
|Tg| . Here, we will merely adopt the final results, which 
are: 

u Q V Q ?/ = sv, n [u n d p hi<f>- 2K al3 u a - 

<S>-\d e - £ N )u ] (21a) 
u„ = Ka^u 13 - u n u a d a In $ + T n , (21b) 
K = h a 0U a u + eu 2 n . (21c) 

where we have defined T a = e a ■ T n = n ■ J 7 , and an 
overdot indicates d/dX. We can express both u a and u n 
in terms of the foliation parameters: 

u a = y a + W a , u n = eM. (22) 

This form of the equations of motion has the virtue of be- 
ing written entirely in terms of tensorial quantities on E^, 
which makes it invariant under n-dimensional coordinate 
transformations. Also note that one of the equations Q21|l 
is redundant; for example, if one contracts (|21a|) with up 
and makes use of (|21c|l . one can recover Q21bfl . 

Now, if a test particle is confined to a given S^ hy- 
persurface, its ^-coordinate must obviously be constant. 
This implies u n = 0, which by (|21b|l yields 

= K af3 u a uP + T n . (23) 

In other words, if the normal force per unit mass 
equals —K a f3U a u 13 , then the particle can be hypersurface- 
confincd. Since this quantity is quadratic in the par- 
ticle's n-velocity, we can identify it as the generalized 
centripetal acceleration in curved space. Indeed, in [l^ 
it was shown that when a particle is confined to the 
world tube of a circle K x S embedded in 3-dimensional 
Minkowski space, T n reduces to the familiar v 2 jr from 
elementary mechanics. 

Now, if one member So of the T,£ foliation satisfies 
K a /3 = 0, then it is obvious that no external centripetal 
force T A is required to ensure confinement. Indeed, when 
the extrinsic curvature vanishes one solution of the freely- 
falling equations of motion is 

y a V a / = 0, 1 = 0; (24) 

i.e., the geodesies of So are also geodesies of M. As 
mentioned in Section [I] surfaces with this property are 
termed totally geodesic and they represent equilibrium 
surfaces for freely- falling test particles. We want to know 
how to tell if such surfaces represent stable or unstable 
equilibria. 

To answer this, we can attempt to linearize the equa- 
tions of motion about So; that is, we consider the motion 
of a test particle very close to the equilibrium hypersur- 
face. To simplify matters, we will adopt the canonical 
gauge l(T§)) discussed above. Then, it is straightforward 
to derive expressions for deh a p and diK a p [jjj: 

dih a p = 2K a/3 , (25a) 
d&Kap = Ka^K^p — E a p, (25b) 



where E a p = e^epn c n D Racbd- This n-tensor can be 
related back to more familiar quantities by suitable ma- 
nipulation of the Gauss-Codazzi relations: 

R AB ei ef = R a /3 + e[E a(} + K a ^{K (}a - Kh Pli )\. (26) 

Now, without loss of generality we can suppose that So 
corresponds to I — 0; from this, it follows that \l\ is 
"small" in our approximation. 3 If we let °/i a /3 denote the 
induced metric on So, °R a f3 the Ricci-tensor on So, etc.; 
then we have 

h a p = O h a[3 + O(f 2 ), (27a) 
K a!3 = e(°R aP - °RABe£ejl)£ + O(f). (27b) 

Furthermore, we suppose that the n-velocity of our par- 
ticle to be approximately tangent to a geodesic on So: 

u a = U a + Su a , U a V a U p = 0. (28) 

Here, Su a is considered to be a small quantity; that is, 
we are really considering perturbations of the confined 
trajectory tangent to U A — U a e A . (We will comment on 
the validity of this assumption below.) Then, to lowest 
order in small quantities, equation (|21b|) reduces to 

e=(°R aP -°R AB e A e^)U a U^+--- (29) 

If S is a stable equilibrium for test particles, we require 
that I/ 1 < 0. This condition translates into the following 
condition for the confinement of test particles on So: 

(°R a p - Q R AB e A e B p )U a UP < 0. (30) 

In order to interpret this result, we define the following 
quantities: 

ee °R AB e A e z U a U?, (31a) 
= °R a0 U a U^. (31b) 

Here, p g n+1 ^ is our definition of the local gravitational 
density of higher-dimensional — or real — matter as 
measured by an observer freely- falling along So- It is 
guaranteed to be positive if the strong energy condition 
is satisfied in the bulk, or at least in the vicinity of the 
totally geodesic surface. Note that it is possible to de- 
fine the gravitational density with a different normaliza- 
tion constant to obtain a "nicer" expression in the per- 
fect fluid case. That is, for a (3 + l)-dimensional perfect 
fluid, an observer comoving with the fluid will measure 



3 Technically, we require that i be small compared to the charac- 
teristic size of the components of the curvature tensors on So 
and M in order to justify dropping the 0(£ 2 ) terms in 1271 . In 
practical terms, this means I should be much less than the radii 
of curvature of both manifolds. 
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Pg 



(4) 



h(p + 3p) using our definition. One might be mo- 



tivated to give a different definition of that does not 
involve the ^ prefactor, but such a choice would unneces- 
sarily complicate the following discussion with spurious 
dimension-dependant terms. 

Now, it is clear that p g n ^ is the lower-dimensional 
cousin of p g n+1 \ but the precise interpretation is a lit- 
tle more subtle. Imagine an observer living on Eo that 
is entirely ignorant of the ^-direction. Assuming that 
this observer can measure the local n-geometry and be- 
lieves in the Einstein field equations, he will interpret the 
Einstein n-tensor of Eo &s being proportional to some 
effective stress-energy tensor. In other words, he will 
conclude that his local n-geometry is determined by an 
effective distribution of matter-energy. From this, it fol- 

(n) 

lows that p g is the gravitational density of the effective 
n-dimensional matter. 

A natural question is: how is the lower-dimensional 
effective matter related to the "real" higher-dimensional 
distribution? It is not hard to derive the following expres- 
sion for the stress-energy tensor on Eo from the Gauss- 
Codazzi relations: 



G a(} = -eE al3 



G AB e^ B 



Gab^ 



AB 



(2=f) Tr(G)] K*, (32) 



where all quantities are understood to be evaluated on 
Eo and we have made note of K a p = 0. The last two 
terms on the right are explicitly related to the higher- 
dimensional stress-energy tensor, but the first can be re- 
written using 

E a p — e A e B n c n D CACBD + n (n-i) £ ^ a l 3 

+ (RA B n A n B h a[ , + sR AB e A e B ) . (33) 

The first term on the right shows how E a p is in part 
determined by the Weyl (n + l)-tensor Cabcd on M. 
This curvature tensor is related to the geometrical or 
gravitational degrees of freedom of the bulk, which are 
not directly fixed by the (n + l)-dimensional field equa- 
tions. This implies that the Einstein n-tensor G a p on AI 
is not entirely determined by the distribution of higher- 
dimensional stress-energy — there is a purely geometric 
contribution from the appropriate projection of Car cd - 
We call this contribution the induced [23] or shadow |2j[ 
matter stress-energy tensor because it represents a source 
of the lower-dimensional Einstein equation that can- 
not be unambiguously attributed to any "real" higher- 
dimensional fields. It follows that the n-dimensional ef- 
fective gravitational density contains contributions from 
both real and induced matter. 

The final point we wish to discuss in this section has 
to do with the validity of our approximations. Recall 
that above, in order to derive equation l|29l) . we assumed 
that 8u a was a small quantity. We now want to describe 
under which circumstances this hypothesis holds by sub- 
stituting the expansion u a — U a + 5u a in equation l|21b[) 



in the canonical gauge with T = 0. For this purpose it 
is useful to assume that U a — U a (y(X)), or equivalently 
dgU a — 0. Coupled with equation l|27a|) . this implies 
that U a S7 'aU 13 — 0(£ 2 ). Under such circumstances, we 
obtain 



-5u fj V g U a , 



(34) 



where we have made use of K a p w — Q E a p x £. It is 
clear that in order to have 8u a be consistently "small", 
we must have that the only term on the right inho- 
mogeneous in 8u a is negligible. Now, since °E a p = 
e(°RABG A e B — ^Rap) we see that its components are of 
the order of the inverse squares of the curvature lengths 
of M and Eo . Then, for our approximations to be valid, 
we need that d(£ 2 )/d\ be small compared to the charac- 
teristic curvature of M or E , whichever is smaller. This 
is a sensible intuitive bound — if either the higher- or 
lower-dimensional manifolds are highly curved we expect 
that true confinement will be difficult to achieve. So, in 
addition to demanding that £ is small in order to justify 
K a f3 ~ —°E a /3 x £, we also need the extra dimensional 
velocity to be relatively tiny. 

To summarize, we have seen that any n-surface Eo 
embedded in M is an equilibrium position for freely- 
falling test particles if it has vanishing extrinsic curva- 
ture. Every geodesic on Eo is automatically a geodesic 
of M, hence the hypersurface is called totally geodesic. 
Furthermore, if a given trajectory confined to E is per- 
turbed off of the equilibrium surface, the acceleration of 
the test particle is towards £ — provided that 

i) the gravitational density of the higher-dimensional 
matter is greater than the gravitational density of 
the effective lower-dimensional matter on E , as 
measured by an observer travelling on the unper- 
turbed trajectory; and, 

ii) both £ and d(£ 2 )/d\ are small compared with the 
characteristic curvature scales of M and Eo. 

Notice that these comments apply to a particular trajec- 
tory only. 



IV. TWO ALTERNATIVE DERIVATIONS OF 
THE STABILITY CONDITION 

The stability condition for freely-falling test particle 
trajectories derived above depended on our decomposi- 
tion of the higher-dimensional equation of motion l|21|) 
and the canonical gauge for the foliation parameters 11911 . 
It is possible to derive the same condition using two dif- 
ferent methods that relax one or both of these assump- 
tions, which is what wc do in this section. 
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A. From the geodesic deviation equation 

Consider a freely-falling test particle on Eo that has 
an (n+ l)-velocity U A = e A U a at a particular instant of 
time. As before, we assume that E is totally geodesic, 
which means that the test particle will remain confined 
on the submanifold for all future times in the absence 
of non-gravitational influences. Now, consider an addi- 
tional test particle that is separated from the first object 
by a vector £ A = £n A . Here, £ is the proper distance 
separating the two particles and is assumed to be small. 
Then, the equation of geodesic deviation says that 

a A = -R A BCD U B fU D , (35) 

where a A is the acceleration of £ A , defined as 

a A = (U ■ V) 2 e A . (36) 

Now, consider n-a. Making use of £ = £n A , we find the 
following expression for this scalar product: 

n-a = e£~ £h AB {U C V c n A ){U D V D n B ). (37) 

Here, we have used gAB = h.AB +£nAns, n A U B V b^-a = 
0, and t = (U ■ V) 2 1 Because U A is parallel to E , the 
second term on the right reduces to —£K a pK ai U ,3 U 1 = 
0. Hence, we have 

£ = -sR A BCDn A U B n c U D £ = ~eE a/ iU a U p £. (38) 

If we now use (|26fl to substitute for E a p, we immediately 
recover our previous result l|29fl. Hence, the test particle 
located just off Eo will be accelerated towards £ = if 
the stability condition from the previous section p g n+1 ^ > 

p^ holds. Notice that we did not assume the canonical 
gauge for this derivation. 

B. From the Raychaudhuri equation 

To establish the stability condition from the Raychaud- 
huri equation, we again consider a freely-falling test par- 
ticle on E with a trajectory 7 tangent to U A = e A U a . 
Consider some small (n — l)-dimensional region 8V n -\ C 
Eo orthogonal to 7 at some given time Ao. We extend 
SV n -i a small distance of £ on cither side of E to de- 
fine an 71-dimensional region SV n . Since U A is tangent 
to Eo, our test particle's trajectory is also orthogonal to 
SV n ■ Now consider a geodesic congruence centered about 
7 and threading every point within 5V n . At the moment 
of interest, we can take each member of the (n + In- 
congruence to be parallel to 7. This means that the 
subset of the total congruence situated on Eo is an n- 
dimensional geodesic congruence on the submanifold. To 
evolve the orthogonal regions forward in time, we imagine 
that each point in SV n is glued to the associated geodesic, 
so the small region deforms in the same manner as the 
congruence. Since the congruence is instantaneously par- 
allel at Aq, we have that ^(Ao) = 0. 



Now in the canonical gauge, the volume of SV n is re- 
lated to the volume of 8V n -i by 

vo\5V n = 2lvo\5V n -x. (39) 

We can define expansion scalars for both the higher- and 
lower-dimensional congruences: 

°n = i ln(voW n ) - V A C/ A , (40a) 
aA 

0„_i = ^ ln(voW n _i) = V a U a . (40b) 
It is easy to see that at A = Ao we have 

d n = £/t+e n - 1 , (41) 

where we have made use of ^(Ao) = 0. The Raychaudhuri 
equation applied to each congruence gives that 

6 n = -V A U B V B U A - RabU a U b , (42a) 
n -i = -V a U p Vf}U a - Q R aP U a U p . (42b) 

A reasonably quick calculation reveals 

V A U B V B U A 
= {h AC + en A n c ){h BD + en B n D )V C U D V B U A 
= V a U p V U a - 2eK aP U l3 e < Xn B V B U A 

+ (U a d a \n<f>) 2 
= V^V^. (43) 

In going from the third to fourth line we have made use 
of the fact that Eo is totally geodesic (K a p — 0) and our 
gauge choice ($ = 1). Putting it all together, we get 

£ = l\p^-pf +1) ). (44) 

The geodesies of the congruence will accelerate towards 
the equilibrium hypersurface if the quantity in the square 
brackets is negative. This yields the same stability con- 
dition as before: pg < 

V. EXAMPLES 

(n) 

We have now established the stability condition p g < 
p^g l+1 ^ using three different methods in general geomet- 
ric backgrounds. In this section, we give a few examples 
of specific manifolds containing totally geodesic hyper- 
surfaces. For each case, we explicitly confirm that our 
stability condition for confined trajectories is correct. 

A. Warped-product spaces 

For our first example of test particle confinement, we 
consider the so-called warped product metric ansatz: 

ds 2 (M) = e n q af} dy a dy 13 + sd£ 2 , (45) 
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where the warp factor f2 = £l(£) is independent of the 
ri-dimensional y coordinates and the warp metric q a p — 
1ap{y) is independent of the extra dimensional £. The 
induced metric on S^ hypersurfaces is h a p = e n q a fj. We 
will assume that the bulk is an Einstein space satisfying 



G 



AB 



AB, 



Rab = 



2A 
i-l 



9AB- 



(46) 



Then, solutions for the warp factor and warp metric are 
easily found [l9| : 



Jl/2 = A 



cos uj£. 



eA > 0, 



coshwi?, eA < 0, 



2AA 2 

= qad, 

n 

.2 2|A| 



n(n 



(47a) 

(47b) 
(47c) 



where A is a constant and R a p is the Ricci tensor formed 
from either the induced or warp metrics (both n-tensors 
give the same result). Essentially, the above states that 
the warp metric can be taken to be any n-dimensional 
solution of the Einstein field equations sourced by a cos- 
mological constant A„ = A 2 A(n — 2)/n. 

Now, for both of the cases shown in (|47ap . it is clear 
that the I — hypersurface S is totally geodesic. 
If we consider some timelike geodesic confined to So, 
the higher- and lower-dimensional gravitational densities 
measured by an observer travelling along that geodesic 



where we have used 



2A 



1' 



rq ! 
J n 



1 = U-U = h aP U a U = A 2 q nf1 U a U f} 



(48) 



(49) 



So, the condition that trajectories on So be stable against 
perturbations in this case simply reads: 



A < 0; 



(50) 



i.e., the bulk must have anti-deSitter characteristics. 4 

This scenario is simple enough to verify our conclusion 
directly from the equations of motion. The particle La- 

r • x, which leads directly 



grangian can be taken as L = ~: 
to the following equation of motion: 



edVL 
2~dl 



(51) 



Now, assuming that our particle is very close to I = 0, 
we can expand d£lj d£ to first order in £ and approximate 



4 It is a misnomer to call the bulk anti-deSitter space in this case, 
it merely satisfies the same Einstein equations as anti-deSitter 
space. 



h a pu a uP 
obtain 

1 



T. In both of the relevant cases eA «s 0, we 



2A 



n(n — 1) 



£ + 0(£ 2 ) = (sgnA)w 2 £ + 0(f 2 ). (52) 



The obvious stability condition from this expression is 
A < 0, which matches the above result precisely. Perhaps 
not surprisingly, the frequency of oscillation about £ = 
is the same frequency found in the warp factor. 

We finish by noting that this metric ansatz easily lends 
itself to toy models of spherically-symmetric astrophys- 
ical situations. For example, suppose that we believed 
that there was a — suitably tiny — cosmological con- 
stant A4 permeating the immediate vicinity of the Sun. 
We could construct a 5- dimensional description by taking 
q a p to be the 4-metric around a A4-black hole: 



q af3 dy a dy f3 

m 



dr 2 

-f(r)dt 2 + ——+r 2 dn 2 , (53a) 



2M 
r 



m 

-A 4 r 2 



(53b) 



where dVl\ is the metric on the unit 2-sphere. Using equa- 



tion 147bl) , we see that A4 is related to the 5-dimcnsional 
cosmological constant by 



A, 



1 Al 



A A. 



(54) 



Interestingly, A4 has the same sign as A. So in this sce- 
nario, test particle like comets, asteroids and planets will 
have a stable equilibrium about £ = if A4 < 0; that is, 
the vacuum energy has negative density. If A4 is positive, 
as suggested by recent cosmological observations, the S 
4-manifold will be gravitationally repulsive to test par- 
ticles. Hence, if a particle at £ = were to acquire a 
small extra dimensional velocity — perhaps by emitting 
gravitational radiation into the bulk — there would be no 
guarantee that it would return to our "native" spacetime. 



B. The Liu-Mashhoon- Wesson metric 

In 0|, test particle trajectories in the following 5- 
geometry were considered: 



ds 2 M) = -b 2 (tj)dt 2 + a 2 {t,£)da 2 k + d£ 2 
a 2 (t,£) = (t 2 + k)£ 2 + K 



b{t,£) 



da 2 



t 2 + k' 
{{t 2 + k) 2 ( 2 - K] 



(t 2 + kf/ 2 [{t 2 + k) 2 £ 2 + K] 1 ' 2 ' 

dr 2 9 , 

1 - 2 dVL 2 2 . 



1 — kr 2 



(55a) 
(55b) 

(55c) 
(55d) 



Here, K, is an integration constant and k = 0, ±1. This 
line element is a special case of one originally discovered 
by Liu & Mashhoon 20] and subsequently re-discovered 
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by Liu & Wesson [21j, and is a solution of the 5- 
dimensional vacuum field equations Rab = 0. It is inter- 
esting because the line element on each of the E^ hyper- 
surfaces is of the cosmological Robertson- Walker form. 
However, this 5-geometry has recently been shown to be 
isometric the 5-dimensional topological Schwarzschild so- 
lution mm. 

It is easy to confirm that £ = is a totally geodesic 
4-surface in the geometry (|55|) with line element: 



^ s (s ) — 



K 



t 2 + k 



dt 2 



{t 2 + kf 



da 2 



(56) 



which can be shown to be isometric to a radiation- 
dominated cosmology. Note that in order to have a sen- 
sible solution, wc need to ensure that 



K 



f 2 + k 



> 



(57) 



by judiciously choosing k and restricting the range of t. 
Now, the tangent vector to a comoving geodesic path on 
S is 



U a d a = (t 2 + k) 



t 2 



-di 



(58) 



The gravitational densities measured by such an observer 
are easily calculated from the basic definitions (1311) : 



P (5) = 



0, 



(4) _ 3(t 2 + fc) 2 



K 



(59) 



Since the bulk is devoid of matter in this case, the stabil- 
ity of the comoving trajectory on So demands p)P < 0. 
In other words, the £ = hypersurface will be gravi- 
tationally attractive only if the strong energy condition 
is violated on So; i- e -, if /C < 0. Notice that in order 
to have a negative value of the integration constant IC, 
the inequality (|57|l implies that k = — 1 and we restrict 

te (-1,1). ' 

This conclusion is odd enough to warrant direct ver- 
ification from the higher-dimensional geodesic equation. 
The 5-dimensional Lagrangian for comoving trajectories 
is 



-b 2 {t,£)i 2 +t 



(60) 



We can obtain an equation for £ by extremizing the ac- 
tion, which yields 



1 -2 d ,2, « ( ^ 

■2 t di b ^ = [wTk 



£ + 0{£ 3 ). 



(61) 



We can use the solution for U a above to approximate 



(r + k) 3 /!C and write 



£ = 



3(t 2 + kf 
K 



£ = P&L 



(62) 



Here, we have omitted the higher-order terms from the 
equation of motion. This is what is expected from Q29JI. 
and confirms to us that the comoving trajectory with 
£ = is stable only if pp < 0. 



VI. SUMMARY AND DISCUSSION: AN 
ENERGY CONDITION FOR HIGHER 
DIMENSIONS 

In this paper, we have shown that confined particle 
trajectories on totally geodesic ro-surfaces embedded in 
(ri + l)-dimensional bulk manifolds are stable against 



(n+l 



00 

> Pg ■ 



Here. 



p g n+1) and 



small perturbations if p g ' 

(n) 

p g are the gravitational densities of M and Eo as mea- 
sured by observers travelling on the confined trajectories, 
respectively. We established this result using a covariant 
decomposition of the higher-dimensional equation of mo- 
tion in Section IlIII the equation of geodesic deviation in 
Section llVAl and the Raychaudhuri equation in Section 
IIVBI In Section [Vl we gave several concrete examples of 
our results involving warped product and Liu-Mashhoon- 
Wesson metrics. 

We conclude by noting that the stability condition as 
formulated above is only applicable to particular geodesic 
paths on So. For some applications, one might want to 
ensure that all the timelike trajectories through some 
region of Eo are stable against perturbations. It is not 
difficult to see how to generalize our previous results to 
satisfy this stronger demand. Consider the following def- 
inition: 

The Confinement Energy Condition: Let So be an 

n-dimensional totally geodesic Lorentzian subman- 
ifold smoothly embedded in (n + l)-dimensional 
bulk M. Also, let °R a /3 be the Ricci-tensor on So 
and °Rab be the Ricci-tensor on M, both evalu- 
ated at a point P E So C M. The confinement 
energy condition at P is 



(°R AB eief-°R a0 )U a U >O, 



(63) 



where U a is an arbitrary timelike vector tangent 
to S . If the confinement energy condition holds in 
some neighbourhood N[P] C So of P, then any test 
particle travelling along a timelike trajectory on So 
passing through P will be stable against small per- 
turbations while in N[P]. 

There are obviously significant similarities between 
this and the familiar strong energy condition from 4- 
dimensional relativity, and we note that it can be used 
to place conditions on the densities and principle pres- 
sures associated with the Einstein tensors of M and 
So- It is clear that for the examples of Section the 
confinement energy condition is satisfied in the warped- 
product Einstein space situation if A < 0, and in the 
Liu-Mashhoon- Wesson metric if the 4-dimensional strong 
energy condition is false on So. We have no doubt that 
it would be interesting to apply this condition to other 
higher-dimensional situations with totally geodesic sub- 
manifolds, but this is the subject for future work. 
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